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Quantum process tomography is a procedure by which an unknown quantum operation can be fully 
experimentally characterized. We reinterpret Choi's proof of the fact that any completely positive 
linear map has a Kraus representation [Lin. Alg. and App., 10, 1975] as a method for quantum process 
tomography. Furthermore, the analysis for obtaining the Kraus operators are particularly simple in this 
. method. 

1 Introduction 

Hi ! 

The formalism of quantum operation can be used to describe a very large class of dynamical evolution of 
quantum systems, including quantum algorithms, quantum channels, noise processes, and measurements. 



> 
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The task to fully characterize an unknown quantum operation £ by applying it to carefully chosen input 
state(s) and analyzing the output is called quantum process tomography. The parameters characterizing the 
quantum operation are contained in the density matrices of the output states, which can be measured using 
quantum state tomography . Recipes for quantum process tomography have been proposed |l2|, |[ [|, [| ||] . 
In earlier methods Jl2|, ||, £ is applied to different input states each of exactly the input dimension of £ . 
In H ||, £ is applied to part of a fixed bipartite entangled state. In other words, the input to £ is entangled 
with a reference system, and the joint output state is analyzed. 

Oh! Quantum processing tomography is an essential tool in reliable quantum information processing, allowing 

error processes and possibly imperfect quantum devices such as gates and channels to be characterized. The 
method in has been experimentally demonstrated and used to benchmark the fidelities of teleportation 
pp| | and the gate CNOT ||, and to investigate the validity of a core assumption in fault tolerant quantum 
computation B. 



The number of parameters characterizing a quantum operation £, and therefore the experimental resources 
for any method of quantum process tomography, are determined by the input and output dimensions of £. 
However, different methods can be more suitable for different physical systems. Furthermore, each method 
defines a procedure to convert the measured output density matrices to a desired representation of £, and a 
simpler procedure will enhance the necessary error analysis. 

In this paper, we describe in detail the method initially reported in H , which is derived as a simple corollary 
of a mathematical proof reported in ||. Our goal is two-fold. We hope to make this interesting proof 
more accessible to the quantum information community, as well as to provide a simple recipe for obtaining 
the Kraus operators of the unknown quantum operation. In the rest of the paper, we review the different 
approaches of quantum operations, describe Choi's proof and the recipe for quantum process tomography in 
Sections ||, ||, and |J. We conclude with some discussion in Section ||. 
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2 Equivalent approaches for quantum operations 



A quantum state is usually described by a density matrix p that is positive semidefinite (p > 0, i.e., all 
eigenvalues are nonnegative) with tr(p) = 1. A quantum operation £ describes the evolution of one state p 
to another p' = £ (p). 

More generally, let Hi and H2 denote the input and output Hilbert spaces of £. A density matrix can be 
regarded as an operator acting on the Hilbert space. [] Let B(TCi) denote the set of all bounded operators 
acting on Hi for i = 1,2. We can consider £(M) for any M G B(Hi) without restricting the domain to 
density matrices. A mapping £ from B(Tti) to B(H2) is a quantum operation if it satisfies the following 
equivalent sets of conditions: 

1. £ is (i) linear, (ii) trace non-increasing (tr(£ (M)) < tr(M)) for all M > 0, and (hi) completely positive. 
The mapping £ is called positive if M > in B(Hi) implies £(M) > in B(H.2)- It is called completely 
positive if, for any auxiliary Hilbert space H a , M > in B(Hi <g> H a ) implies (£ ®I){M) > in 
B{H.2 ® H a ) where l a is the identity operation on B(H a ). 

2. £ has a Kraus representation or an operator sum representation (l^, |3[ 

£{M) = Y J A k MA\ (1) 

k 

where Ylk-^k-^-k < I, and / is the identity operator in B(Hi). The Ak operators are called the Kraus 
operators or the operation elements of £. 

3. £{M) = Tr Q [ U(M ® p a )U^ (I ® P Q ) ] . Here, p a € S(W Q ) is a density matrix of the initial state of the 
ancilla, I £ B{H.2) is the identity operator, H2 ®H = Hi ®H a , P G B(H a ) is a projector, and Tr Q is 
a partial tracing over 7i D . 

Each set of conditions represents an approach to quantum operation when the input is a density matrix 
(M = p) . The first approach puts down three axioms any quantum operation should satisfy. The completely 
positive requirement states that if the input is entangled with some other system (described by the Hilbert 
space H a ), the output after £ acts on Hi should still be a valid state. The third approach describes system- 
ancilla (or environment) interaction. Each of these evolutions results from a unitary interaction of the system 
with a fixed ancilla state p a , followed by a measurement on a subsystem H with measurement operators 
{P ,I — P }, post-selection of the first outcome, and removing H Q - The fact that the third approach 
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1 Even though the density matrix can be viewed as an operator, it is important to remember that it represents a state and 
thus evolves as a state rather than an operator. For this reason we use the name density matrix rather than the density operator. 
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is equivalence to the first is nontrivial - the evolutions described by the third approach are actually all 
the mappings that satisfy the three basic axioms. Finally, the second approach provides a convenient 
representation useful in quantum information theory, particularly in quantum error correction (see jfl| for 
a review) . 

Proofs of the equivalence of the three approaches are summarized in ^ There are four major steps, 
showing that the 1st set of conditions implies the 2nd set and vice versa, and similarly for the 2nd and 3rd 
sets of conditions. The most nontrivial step is to show that every linear and completely positive map has a 
Kraus representation, and a proof due to Choi [Bj for the finite dimensional case will be described next. 



3 Choi's proof 



The precise statement to be proved is that, if £ is a completely positive linear map from B(Hi) to 6(7^2), then 

£(M) = J2k A-kMA^ for some 77,2 xni matrices A^, where rij is the dimension of Hi- Let |$) = -7= ^\ 

be a maximally entangled state in Hi ®Hi- Here, {\i)}i=i,..., ni is a basis for Hi. Consider (X(g>£)(M) where 



(2) 



M is an n\ x m array of n± x n\ matrices. The block is exactly 



M = 
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(3) 
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When I® £ is applied to M, the block becomes £(\i)(j\), and 



(T<g>£)(M) = 



/ 1 


\ 





1 






/ 1 


\ 




1 





















(4) 



which is an n\ x ri\ array of n 2 X n 2 matrices. 

We now express (T (g) £){M) in a manner completely independent of Eq. (||). Since M is positive and £ is 
completely positive, (1®£){M) is positive, and can be expressed as {X®£ ){M) — J2?li 2 l a fe}(afc|, where \ak) 
are the eigenvectors of (T(g£)(M), normalized to the respective eigenvalues. One can represent each \ak) as 
a column vector and each (afc| as a row vector. We can divide the column vector \ak) into ri\ segments each 
of length ri2, and define a matrix with the i-th column being the i-th segment, so that the i-th segment 
is exactly Ak\i). Then 



(I®£)(M) = Ek 



(ii4 
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and 



A k \l)(l\A{ 


A k \l)(2\A\ 




A k \l){m\A\ 


A k \2)(l\A\ 


A k \2){2\A\ 




A h \2){ni\A\ 










A k \ ni )(l\A\ 


A k \nx){2\A\ 




A k \ni){m\A\ 



(5) 



Comparing Eqs. (|) and (|) block by block £{M) = J2k A k MA\ for MM = \i)(j\, and thus MM G B(Hi) by 
linearity. 



4 Recipe for quantum process tomography 

The basic assumptions in quantum process tomography are as follows. The unknown quantum operation, £, 
is given as an "oracle" or a "blackbox" one can call without knowing its internal mechanism. One prepares 
certain input states and measures the corresponding output density matrices to learn about £ systematically. 
The task to measure the density matrix of a quantum system is called quantum state tomography jlj . To 
obtain a Kraus representation for £, one needs an experimental procedure that specifies the input states to 
be prepared, and a numerical method for obtaining the Kraus operators from the measured output density 
matrices. 

A method follows immediately from the proof in Sec. [|. We retain all the previously defined notations. The 
crucial observation is that and -^(I®£)(M) correspond to the input and output physical states |$)(<I>| 
and (1 (£> £)(|$)($|) which can be prepared and measured. The procedure is therefore to: 

1. Prepare a maximally entangled state |<I>) in Tii <g> T~L\- 

2. Subject one system to the action of £, while making sure that the other system does not evolve. 

3. Measure the joint output density matrix (2 <g> £ )(|$}($|) = ^-{Z <8> £){M), multiply by m, obtain the 
eigen-decomposition J^ fc \ak)(ak\- Divide \a k } (of length nxfio) into ri\ equal segments each of length 
rii. Ak is the ni x n\ matrix having the i-th segment as its i-th column. 

The maximally entangled state in the above procedure can be replaced by any pure state with maximum 
Schmidt number, \<f>) — a.i(U\i))®{V\i)) where > are real and af — 1. The output density matrix 
p ou t is equal to (I ® £)(\<j>) (<fi\) = J2 iyj aiOLj(U\i) (j\U^) ® £{V\i){j\V^). One computes (l/t ®/)p out (tT®i), 
divides the block by aijay, and performs eigen-decomposition to obtain a set of A k operators. The 

Kraus operators of £ are given by A k V'. 
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5 Discussion 



We have provided an experimental and analytic procedure for obtaining a set of Kraus operators A k for an 
unknown quantum operation. The set of Ak is called "canonical" in ||, meaning that the Ak are linearly 
independent. We remark that any other Kraus representation can be obtained from Ak using the fact that 
£{P) = Ek^kpAi = EkBkpBl if and only if A k = V. u k jB k when u k j are the entries of an isometry 
||. Alternatively, one can replace the eigen-decomposition of (T ® by any decomposition into a 

positive sum to obtain other valid sets of Kraus operators. 

Previous methods of quantum process tomography |l^, [|, |B| involve preparing a set of physical input states 
Pi that form a basis of B(TL\), and measuring £{pi) to determine £. The input states pi are physical 
states, and cannot be chosen to be trace orthonormal, causing complications in the analysis. In contrast, 
the output state in the current method automatically contains complete information on £ for the 
unphysical orthonormal basis (see Eq. (|j)), which greatly simplifies the analysis to obtain the Kraus 

operators. However, the current method requires the preparation of a maximally entangled state and the 
ability to stop the evolution of the reference system while £ is being applied. The previous methods are 
more suitable in implementations such as solution NMR systems, while the current method is more suitable 
for implementations such as optical systems. 

Any efficient quantum process tomography procedure consumes approximately the same amount of resources, 
which is determined by the number of degrees of freedom in the quantum operation. In general, to measure 
an n x n density matrix, n 2 ensemble measurements are needed, requiring w 0{n 2 ) steps. The previous 
methods require the determination of n\ density matrices each n% x ni and take w 0(n\n\) steps. The 
current method requires the determination of one n\ni x n\n2 density matrix which also requires ~ 0(nfn§) 
steps. In both cases, the number of steps is of the same order as the number of degrees of freedom in the 
quantum operation and are optimal in some sense. 
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